*DECK USER
      SUBROUTINE USER
C================================================================USER
C     THIS ROUTINE PROVIDES USER SPECIFIED GRID INFORMATIONS
C----------
C ----- could be changed into input file Y.P. WAN !!!
C     CALLED BY SETUP
C----------
C     OUTPUT: X,Y,Z,SGSL
C================================================================USER
      INCLUDE 'COML.h'
C
C================================================================USER
      ENTRY USERX
C----------------------------------------------------------------USER
C     THIS ROUTINE SHOULD PROVIDE X(I) VALUES AS USER SPECIFIED.
C----------------------------------------------------------------USER
c----- METCO 9MB, FGM, Pseudo 3-D Particle
      x( 1)= 0.00000D+00
      x( 2)= 3.90000D-02
      x( 3)= 7.80000D-02
      x( 4)= 1.17000D-01
      x( 5)= 1.56000D-01
      x( 6)= 1.95000D-01
      x( 7)= 2.34000D-01
      x( 8)= 2.73000D-01
      x( 9)= 3.12000D-01
      x(10)= 3.51000D-01
      x(11)= 3.90000D-01
c---- METCO 9MB nozzle radius = x(11)
      x(12)= 4.29000D-01
      x(13)= 4.69672D-01
      x(14)= 5.12088D-01
      x(15)= 5.56322D-01
      x(16)= 6.02453D-01
      x(17)= 6.50562D-01
      x(18)= 7.00733D-01
      x(19)= 7.53055D-01
      x(20)= 8.07621D-01
      x(21)= 8.64525D-01
      x(22)= 9.23870D-01
      x(23)= 9.85759D-01
      x(24)= 1.05030D+00
      x(25)= 1.11761D+00
      x(26)= 1.18781D+00
      x(27)= 1.26101D+00
      x(28)= 1.33735D+00
      x(29)= 1.41697D+00
      x(30)= 1.50000D+00
c---- METCO 9MB Torch face outside radius = x(30)
      x(31)= 1.58303D+00
      x(32)= 1.67016D+00
      x(33)= 1.76160D+00
      x(34)= 1.85756D+00
      x(35)= 1.95827D+00
      x(36)= 2.06395D+00
      x(37)= 2.17485D+00
      x(38)= 2.29124D+00
      x(39)= 2.41338D+00
      x(40)= 2.54155D+00
      x(41)= 2.67606D+00
      x(42)= 2.81722D+00
      x(43)= 2.96535D+00
      x(44)= 3.12081D+00
      x(45)= 3.28395D+00
      x(46)= 3.45516D+00
      x(47)= 3.63482D+00
      x(48)= 3.82337D+00
      x(49)= 4.02123D+00
      x(50)= 4.22888D+00
      x(51)= 4.44679D+00
      x(52)= 4.67547D+00
      x(53)= 4.91545D+00
      x(54)= 5.16729D+00
      x(55)= 5.43158D+00
      x(56)= 5.70894D+00
      x(57)= 6.00000D+00
C----------------------------------------------------------------USER
      X(0)=X(1)-X(2)+X(1)
      X(NXT)=X(NXT-1)+X(NXT-1)-X(NXT-2)
      RETURN
C================================================================USER
      ENTRY USERY
C----------------------------------------------------------------USER
C     THIS ROUTINE SHOULD PROVIDE Y(I) VALUES AS USER SPECIFIED.
C----------------------------------------------------------------USER
c----- METCO 9MB, FGM, Pseudo 3-D Particle
      y( 1)=  .00000D+00
      y( 2)= 4.00000D-02
      y( 3)= 8.18216D-02
      y( 4)= 1.25548D-01
      y( 5)= 1.71265D-01
      y( 6)= 2.19065D-01
      y( 7)= 2.69041D-01
      y( 8)= 3.21293D-01
      y( 9)= 3.75924D-01
      y(10)= 4.33044D-01
      y(11)= 4.92765D-01
      y(12)= 5.55205D-01
      y(13)= 6.20489D-01
      y(14)= 6.88746D-01
      y(15)= 7.60111D-01
      y(16)= 8.34726D-01
      y(17)= 9.12740D-01
      y(18)= 9.94306D-01
      y(19)= 1.07959D+00
      y(20)= 1.16875D+00
      y(21)= 1.26197D+00
      y(22)= 1.35944D+00
      y(23)= 1.46135D+00
      y(24)= 1.56790D+00
      y(25)= 1.67931D+00
      y(26)= 1.79578D+00
      y(27)= 1.91756D+00
      y(28)= 2.04489D+00
      y(29)= 2.17801D+00
      y(30)= 2.31719D+00
      y(31)= 2.46272D+00
      y(32)= 2.61487D+00
      y(33)= 2.77395D+00
      y(34)= 2.94028D+00
      y(35)= 3.11418D+00
      y(36)= 3.29600D+00
      y(37)= 3.48609D+00
      y(38)= 3.68485D+00
      y(39)= 3.89266D+00
      y(40)= 4.10993D+00
      y(41)= 4.33709D+00
      y(42)= 4.57460D+00
      y(43)= 4.82293D+00
      y(44)= 5.08257D+00
      y(45)= 5.35402D+00
      y(46)= 5.63784D+00
      y(47)= 5.93459D+00
      y(48)= 6.24485D+00
      y(49)= 6.56924D+00
      y(50)= 6.90840D+00
      y(51)= 7.26301D+00
      y(52)= 7.63376D+00
      y(53)= 8.02140D+00
      y(54)= 8.42670D+00
      y(55)= 8.85045D+00
      y(56)= 9.29349D+00
      y(57)= 9.75672D+00
      y(58)= 1.02410D+01
      y(59)= 1.07474D+01
      y(60)= 1.12768D+01
      y(61)= 1.18304D+01
      y(62)= 1.24091D+01
      y(63)= 1.30143D+01
      y(64)= 1.36469D+01
      y(65)= 1.43084D+01
      y(66)= 1.50000D+01
C----------------------------------------------------------------USER
      Y(0)=Y(1)-Y(2)+Y(1)
      Y(NYT)=Y(NYT-1)+Y(NYT-1)-Y(NYT-2)
      RETURN
C================================================================USER
      ENTRY USERZ
C----------------------------------------------------------------USER
C     THIS ROUTINE SHOULD PROVIDE Z(I) VALUES AS USER SPECIFIED.
C----------------------------------------------------------------USER
      Z(0)=Z(1)-Z(2)+Z(1)
      Z(NZT)=Z(NZT-1)+Z(NZT-1)-Z(NZT-2)
      RETURN
C================================================================USER
      ENTRY USERSL
C----------------------------------------------------------------USER
C     THIS ROUTINE SHOULD PROVIDE SUB-GRID LENGTH SCALE SGSL
C     WHEN ITURB=1
C------------------------------ SGSL=2*MINIMUM DELTA X,Y,Z ------USER
C     WHEN THERE ARE CELLS WITH BAD ASPECT RATIO, THE SMALLEST SIDE
C     WILL BE CHOSEN. THIS IS NOT A GENERAL FEATURE, SO USER SHOULD
C     EXAMINE APPROPRIATE ONE FOR HIS OWN PROBLEM.
C----------------------------------------------------------------USER
      DO 600 IC=1,NXYZT
      IF(NX.EQ.1) THEN
        SGSL(IC)=TWO*MIN(DY(IC),DZ(IC))
        IF(NY.EQ.1) SGSL(IC)=TWO*DZ(IC)
        IF(NZ.EQ.1) SGSL(IC)=TWO*DY(IC)
      ELSEIF(NY.EQ.1) THEN
        SGSL(IC)=TWO*MIN(DX(IC),DZ(IC))
        IF(NZ.EQ.1) SGSL(IC)=TWO*DX(IC)
      ELSEIF(NZ.EQ.1) THEN
        SGSL(IC)=TWO*MIN(DX(IC),DY(IC))
      ELSE
        SGSL(IC)=TWO*MIN(DX(IC),DY(IC),DZ(IC))
      END IF
       sgsl(ic)=3.75d0*sqrt(dx(ic)*dy(ic))
  600 CONTINUE
C================================================================USER
      RETURN
      END
